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An analytical model of dispersive pulse propagation in semi-infinite cylinders due to transient 
axially symmetric end conditions has been experimentally investigated. Specifically, the dispersive 
propagation of the first axially symmetric longitudinal mode in thin wire waveguides, which have 
ends in butt contact with longitudinal piezoelectric ultrasonic transducers, is examined. The method 
allows for prediction of a propagated waveform given a measured source waveform, together with 
the material properties of the cylinder. Alternatively, the source waveform can be extracted from 
measurement of the propagated waveform. The material properties required for implementation of 
the pulse propagation model are determined using guided wave phase velocity measurements. Hard 
tempered aluminum 1100 and 304 stainless steel wires, with 127, 305, and 406 fim diam, were 
examined. In general, the drawn wires were found to behave as transversely isotropic media. 

PACS numbers: 43.35. Pt, 43.20.Jr 


INTRODUCTION 

Historically, pulse propagation in thin cylindrical 
waveguides has been of interest for applications including 
the material testing of wires and fibers, 13 ultrasonic delay 
lines, 4 medical ultrasonics, 5 and embedded sensor materials. 6 
Propagation of ultrasonic waves in cylinders has been the 
subject of several investigations, both analytical and 
experimental. 7 ” 20 In high-frequency and broadband applica- 
tions, these cylinders may exhibit considerable dispersion 
when an acoustic pulse propagates through the waveguide. 
Here, dispersion is the dependence of acoustic guided wave 
phase velocity on frequency, resulting in distortion of the 
waveform of a transient plane-wave pulse depending on the 
axial location. This pulse distortion can make the interpreta- 
tion of measured waveforms difficult. The present paper 
combines experimental measurements and analytical model- 
ing to yield a quantitative and integrated approach to address 
the problem of pulse dispersion in cylindrical waveguides. 
The method presented also has significant relevance to the 
problem of material characterization of these waveguides. 

I. PULSE PROPAGATION 

While exact 7 0 and approximate 10,11 solutions have been 
used to study the axially symmetric guided wave phase ve- 
locity versus frequency relationship in solid cylindrical 
waveguides, and experimental dispersion measurements 
have corroborated the theory, 12 ' 15 few studies have been 
done to assess the capability these solutions yield for mod- 
eling actual pulse propagation. The earliest analytical solu- 
tion, using exact theory, to the problem of axially symmetric 
pulse propagation in a cylindrical waveguide due to applica- 
tion of a transient end condition can be attributed to 
Skalak. 16 A special solution was obtained to the problem of a 
solid isotropic cylinder with a velocity imparted uniformly to 


the cylinder end in the form of a step function. Subsequently, 
Folk et aO 1 used exact theory to obtain a solution for a solid 
isotropic cylinder with arbitrary time harmonic mixed (com- 
bined displacement and stress) end boundary conditions. 
Moreover, an attempt was made by Folk et al. ]1 to experi- 
mentally validate the solution in terms of its ability to model 
pulse propagation. Studies by Miklowitz and Nisewanger, 18 
Hseih and Kolsky, 19 and McNiven and Mengi 20 have at- 
tempted to experimentally validate propagation models 
based on approximate theories. Only qualitative agreement 
between theoretical and measured waveforms was obtained 
in these studies. Furthermore, the pulses examined in all of 
these studies were dominated by lower-frequency propaga- 
tion of the first longitudinal mode. The pulses, resulting from 
either ideal or modified (some finite rise time) step loading, 
essentially consist of a low-frequency leading head and a 
higher-frequency dispersive tail. While the current study also 
addresses first longitudinal mode propagation, the transient 
end conditions applied give rise to more varied dispersive 
pulse distortion. The frequency range considered extends 
from low frequency, where the fundamental tends toward 
nondispersiveness, to approaching the second longitudinal 
mode cutoff, where the fundamental becomes highly disper- 
sive. An integrated analytical and experimental study of iso- 
tropic hollow cylinders using the same numerical method 
and similar measurement techniques 21 involves the propaga- 
tion of higher-order longitudinal modes. 


II. PROPERTY MEASUREMENT 

Because the dispersion of the guided waves depends on 
the elastic properties, cylindrical guided wave propagation 
has been used for the material characterization and testing of 
wires and fibers. Meeker and Meitzler^ suggested several 
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methods for determining material properties by measuring 
acoustic guided wave velocities, and a thorough review of 
the subject has been given by Curtis. 1 

Principally, velocity measurements of the first torsional 
mode and the long wavelength limit of the first longitudinal 
mode have been used to obtain the shear and Young’s 
moduli, respectively, for isotropic cylinders. 1,2 The first axi- 
symmetric torsional mode is nondispersive, propagating at 
the shear wave velocity (G7p) 1/2 . In the long wavelength 
(low-frequency) limit, the velocity of the first longitudinal 
axisymmetric mode tends toward the bar velocity (E/p) lu . 
As the limit is approached, group and phase velocities tend 
toward a constant and equal value, indicating nondispersive 
behavior. 

However, the drawing process in the manufacture of 
wires and fibers can affect the material microstructure, fre- 
quently causing the material to be anisotropic. 1 Transverse 
isotropy, with five independent elastic constants, is consid- 
ered to better represent the state of the material than isotropy. 
The five independent elastic properties for transversely iso- 
tropic media are taken as £ LL , G lx , where 

E, G , and v are Young’s moduli, shear moduli, and poisson 
ratio, respectively, L is longitudinal, and T is transverse. 
Moreover, the propagation of the fundamental longitudinal 
mode, within the frequency range considered, does not de- 
pend on the transversely isotropic property to which the low- 
est torsional mode is sensitive, G LX . Rather, the strongest 
dependence is on E LL and y LX . 

Longitudinal mode measurement to obtain £ LL can be 
problematic. When the velocity of the long wavelength limit 
is used, a transient pulse with sufficient low-frequency con- 
tent is propagated and its first arrival measured. Since the 
longest wavelengths of the first longitudinal mode propagate 
with greatest velocity, this time-of-flight measurement yields 
the velocity of the long wavelength limit. However, the lead- 
ing head of a dispersive pulse often rises very gradually mak- 
ing accurate assessment of when the first arrival occurs dif- 
ficult, as discussed by Curtis, 1 and evident in the low- 
frequency experimental waveforms of other studies. 17-20 
Consequently, the measurement may be subject to some er- 
ror. Instead, the time shift between the peak of the wave 
envelope of a propagated broadband pulse is often matched 
with that of the initial pulse. 2 However, the frequency and 
wavelength of such a velocity measurement is not at all 
clear. 24 

Alternatively, phase or group velocity of particular cy- 
lindrical guided wave modes can be measured over a range 
of frequency. Material properties are then varied to obtain 
correspondence between theoretical velocities and the disper- 
sion measurements. As noted by Meeker and Meitzler, 22 dis- 
persion measurements of the first axially symmetric longitu- 
dinal mode are sufficient to find both elastic constants for 
isotropic media. As will be illustrated here, combination of 
such measurement with measures of first torsional mode ve- 
locity yields £ LL , v LX , and G LX for transversely isotropic 
media. The extent of the material anisotropy can then be 
gauged. 

Phase velocity measurement methods are preferred to 
group velocity measurement methods since the wave propa- 
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FIG. 1. Experimental apparatus. 

gation is generally dispersive. Inherent difficulty in group 
velocity measurements for dispersive media has been 
noted. 24,25 The problem stems from the distortion of the 
propagated waveform relative to the initially generated 
waveform and the resulting ambiguity in correlating the two 
different waveforms, and also, in determining the frequency 
for which the group velocity is being measured. Correspond- 
ingly, precision and accuracy for phase velocity measure- 
ment can be expected to exceed that obtained for group ve- 
locity measurement for dispersive pulse propagation. The 
scatter of group velocity measurements which commonly oc- 
curs for propagation in dispersive media can be seen in a 
recent study of sound in wires by Nicholson et a/. 15 The 
precision of the measurement and smoothness of the curve 
determined by these measurements is critical to accurate and 
repeatable material property characterization. The elastic 
properties determined from the phase velocity measurements 
are used here as input to the pulse propagation model. 

III. WAVE EXCITATION AND MEASUREMENT 

The method by which elastic pulses are launched into 
the wire specimens and measured after propagating the 
length of the sample is illustrated in Fig. 1. Similar tech- 
niques have been used in other studies of ultrasound propa- 
gation in wires. 15 As shown in Fig. 1, matched pairs of pi- 
ezoelectric longitudinal or shear wave transducers are placed 
in butt contact with each end of the wire and operate in pitch 
catch mode. For this study, 5 MHz center frequency pana- 
metrics V109 (longitudinal) and V156 (shear) transducers 
were used. Glycerin was used as couplant to aid in sound 
transmission across the transducer/wire interface. The source 
or excitation transducer was provided a gated sine wave or 
toneburst electrical signal. Because of the thinness of the 
wires, the efficiency at which sound could be introduced into 
the waveguides was very low. Consequently, gain in the re- 
ceiving system between 60 and 100 dB was required and 
signal averaging between 300 and 600 data records was per- 
formed to improve signal to noise. Signals were captured and 
averaged on a digital oscilloscope and processed on a PC. 

With longitudinal wave transducers in use, the pulse 
launched in the medium would be composed primarily of the 
first axially symmetric longitudinal mode. However, depend- 
ing on factors such as the wire thickness, excitation fre- 
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FIG. 2. Measured propagated pulse in a 24.76-cm-long, 12-mil-diam alumi- FIG. 3. Measured propagated pulse in a 25.22-cm-long, 5-mil-diam 304 

num IKK) wire due to excitation of longitudinal piezoelectric source trans- stainless steel wire due of excitation to shear piezoelectric source transducer 

duccr with l -MHz (top) and 5-MHz (bottom) tonebursts. with 700-kHz (top) and 1.5-MHz (bottom) tonebursts. 


quency, couplant effects, and obviously the contact between 
source transducer and wire, the first-order flexural mode was 
also excited. In general, the flexural mode propagates much 
more slowly than the axisymmetric longitudinal mode. Con- 
sequently, the waveform due to each of the two modes can be 
temporally separated after some length of propagation. How- 
ever, with increasing frequency, the dispersion of the wave- 
guide is such that the first longitudinal mode is observed to 
propagate more slowly while the flexural mode is observed 
to propagate faster. Propagation of higher-frequency harmon- 
ics may result in overlap and interference between axially 
symmetric longitudinal and nonsymmetric flexural waveform 
components. The longitudinal and flexural components of 
typical measured waveforms for bursts with increasing cen- 
ter frequency are shown in Fig. 2. Longitudinal wave trans- 
ducers are used for generating and receiving the wave. The 
waveform component associated with the flexural mode was 
eliminated through windowing for both the phase velocity 
measurements and the extraction of the source signal from a 
measured received waveform. 

Torsional modes were excited by using shear wave 
transducers and intentionally placing only an edge of the 
wire end into contact with the source transducers face. As a 
result of such contact, both longitudinal and torsional axi- 


symmetric modes and at least first-order flexural modes were 
excited. By varying the frequency of the toneburst excitation 
of the shear wave transducer, and observing which wave 
packet remained fixed with regard to arrival and did not dis- 
tort relative to the source waveform (nondispersive), the tor- 
sional wave packet could be identified. Torsional wave pack- 
ets are shown for bursts of 0.7 and 1.5 MHz in Fig. 3. Since 
dispersive first-order flexural waves are also excited, caution 
must be used as flexural waves will often propagate at simi- 
lar velocities relative to torsional waves. Consequently there 
is overlap of the waveforms over broad ranges of frequency. 
In Fig. 3, the faster propagation of the flexural mode when a 
higher-frequency burst was used, results in the beginning of 
interference between the nondispersive torsional mode and 
the highly dispersive flexural mode. 

Source measurements were performed by placing the 
matched transducers into face to face contact without the 
intervening wire. It should be noted that source measurement 
by this method is imperfect since it neglects all transducer/ 
wire interface effects. In particular, some variation in the 
behavior of the source transducer when electrically excited in 
unloaded (with wire contacting transducer face) as opposed 
to loaded (matched transducer face-to-face contact) condition 
could be anticipated. The effect of multiple reflections be- 
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tween the two transducer faces on the source signal measure- 
ment was also neglected. 

IV. BOUNDARY CONDITIONS AND ANALYTICAL 
MODELING 

The cylinder end boundary conditions used in the ana- 
lytical modeled, implemented here, approximate the test con- 
figuration of Fig. 1. The longitudinal source transducer is 
modeling as a piston source in smooth contact with the cyl- 
inder end. The boundary conditions applied are a time de- 
pendent pure axial pressure pulse p(t) and vanishing shear 
stresses acting on the cylinder end surface. On the surface at 
z — 0 then, normal traction is uniform [p z = p(t)] and shear 
tractions are zero (p T ~p^— 0). Because the entire cylinder 
end is idealized to be in uniform contact with the longitudi- 
nal transducer, only axially symmetric mode excitation and 
propagation is modeled. Since flexural waves are excited ex- 
perimentally, as was shown in Fig. 2, the model boundary 
conditions do not fully characterize the actual physical end 
conditions. However, assuming the end conditions can be 
decomposed into axially symmetric and nonsymmetric parts, 
it should be of great utility to develop a capability for both 
receiver waveform prediction and source waveform extrac- 
tion based on the axially symmetric wave propagation, with- 
out necessarily characterizing the nonsymmetric part of the 
boundary conditions. 

In a previous paper, 26 a solution for the problem of a 
time varying uniform (across the end face) axial displace- 
ment pulse in combination with (1) vanishing shear stresses 
and (2) vanishing circumferential and radial displacements, 
was given. The analytical development given in Kohl 
et al 2b is generalized in the current paper to boundary con- 
dition problems involving any well posed combination of 
prescribed axially symmetric, transient distributed forces or 
discretized displacements, specified at the cylinder end. The 
specific aforementioned boundary conditions, modeling the 
butt contact of a longitudinal ultrasonic transducer, give rise 
to a pure stress boundary condition. 

The numerical method to be described has diverse ap- 
plication. The method is readily applicable to problems in- 
volving both solid and hollow cylindrical waveguides. Fur- 
thermore, the media of the waveguides need not be isotropic 
or homogeneous. A solution for problems involving aniso- 
tropic and radially inhomogeneous cylinders, such as fila- 
ment wound composites, can be obtained. For conciseness, 
the development here follows directly from that given in 
Kohl et a /., 26 with only new formulation presented. Discus- 
sion in Kohl et ai 2b also provides further details on imple- 
mentation of the modeling. 

The eigenproblem of Eq. (3) in Kohl et al 2b constitutes 
numerical solution for both the mode phase velocity disper- 
sion relation, and the mode shapes or particle displacement 
distribution for each mode within the plane of the advancing 
wavefront. The integrations involved in the formation of the 
matrices of Eq. (3) 26 can be evaluated exactly except for [t 3 ] 
which has a singularity at the center of a solid cylinder. Ei- 
ther Gaussian numerical integration or introduction of a pin- 
hole to the center of the cylinder are satisfactory techniques 
for dealing with this singularity. 


A variational formulation is used to determine the end 
condition equation governing the modal amplitudes A n . Ap- 
plication of the principle of virtual work yields 

J ( z Su r + H+ (T zz Su z )dA 



( p r 8u r +p e 8u e +p z 8u z )dA , 


( 1 ) 


where p t are the prescribed distributed forces acting on the 
cylinder end, u x are the displacements, and a l} are the 
stresses. Equation (1) is applied to each element in the radi- 
ally discretized system, yielding for the kih element 

Jr b 


Xdr{A} = [SA] T [Qf k ' T f rf [N] T {p}r dr, (2a) 

Jr b 


where 


[C*] {k) = 


0 0 


14 

C 2 4 

13 

^23 


0 

0 

C 55 

C 56 ' 

k 

^34 

C44 

0 

0 

(2b) 

C 33 

C 34 

0 

0 _ 



and the distributed force vector {p} is 


[pY=iPrPePz\- 


(2c) 


For definition and explanation of the quantities appearing in 
Eqs. (2a)-(b) the reader is again referred to Kohl et ai 2ti The 
partition of the constitutive matrix in Eq. (2b), corresponding 
to the constitutive relations of Eq. (9) in Kohl et al., 26 are for 
“cylindrically anisotropic ” 27 media. This is the most general 
form of the constitutive relations considered. For materials 
which are transversely isotropic in the natural coordinate sys- 
tem of the cylinder, as the wires are considered here, 
^56 = ^34~^24 = Ci 4“0 and C 13 = C 2 3 . 

For a pure axial pressure pulse, the impulse response is 
obtained using the boundary condition equation 

{p} T ={ 001 ]. (3) 


The integral on the left side of Eq.(2a) is the element 
contribution to an amplitude stiffness matrix. The integral on 
the right side of Eq. (2a) is the element contribution to a 
consistent nodal force vector. Assembly of the global stiff- 
ness matrix and consistent nodal force vector are done in the 
usual manner. The integrations of Eq. (2a) can be performed 
exactly for both solid and hollow cylinders. 

As with traditional finite elements, either force or dis- 
placement may be specified on a node-by-node basis for each 
degree of freedom. The consistent nodal force vector repre- 
sents a variational argument for node-by-node specification 
of force given prescribed surface traction. A solution using 
an unaltered amplitude stiffness matrix and consistent nodal 
force vector addresses the problem of a pure stress boundary 
condition problem. Such a solution is being examined here. 
Alternatively, using Eq. (2), in Kohl et al. 26 evaluated at 
z— 0, the nodal displacements for the kth element are given 
by 
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(6a) 


(4) 

Nodal displacement is prescribed by replacing the row of the 
assembled global amplitude stiffness matrix corresponding to 
the displacement degree of freedom being specified with the 
appropriate row of [Q] {k) and replacing the same row entry 
of the consistent nodal force vector with the specified dis- 
placement. 

For specific <o , imposition of the boundary conditions of 
Eq. (3) and solution for the complex mode amplitudes 
[A(a>)] yields the impulse response G(a j) to a pure pressure 
pulse for any of the displacements, strains or stresses of in- 
terest. The time domain responses for any field variable, say 
x(/), may be written as a convolution between the time do- 
main impulse response g(t) and the excitation that here is 
the prescribed pressure p(t) 

x(t)= I p(\)g(t~k)d\. (5) 

Jo 

The frequency domain transform of the field variable x(/) is 
then given as X{ oj) = G(co)P(co). 

Considering the experimental configuration of Fig. 1, the 
analytical model is used to predict the waveform of the axial 
force resultant [/ z = *(/)], or integral of the normal stress or zz 
at the specimen end z = L, for comparison with experimental 
measurement. Given waveguide geometry and properties, the 
impulse response [G(oj)] for f z is computed at discrete fre- 
quencies over sufficient range so that all products 
G(w)P(oj) with significant amplitude are included in the 
computation. As discussed, the source waveform p(t) is ob- 
tained by holding the generating and receiving transducers in 
face-to-face contact without the intervening waveguide. Fast 
Fourier transform (FFT) techniques are used to obtain P(a j). 
An inverse FFT of the product F 2 (tu) = G(to)P(w) yields 
f z = x (t). Or the source p(t) can be extracted from the mea- 
sured axial force resultant f z using an inverse FFT of the 
quotient P(co) = F z (co)/G((o). 

V. PHASE VELOCITY MEASUREMENT AND 
MATERIAL PROPERTY CHARACTERIZATION 

Effective implementation of the pulse propagation 
model requires accurate material characterization. The 
method of phase spectral analysis, given by Sachse and 
Pao, 25 is used with slight modification to measure phase ve- 
locity of the first longitudinal axisymmetric mode in thin 
wires over a 1-6 MHz range of frequency. Together with 
measurement of the first axially symmetric torsional mode, 
three of the five independent material properties for trans- 
versely isotropic media could be determined. 

Measurement of the cylindrical guided wave velocity by 
phase spectral analysis requires measurement of an ultra- 
sonic field at two axial locations. Assuming a single mode is 
associated with the measured signal, an ultrasonic field, 
x(/,z), due to plane harmonic waves, can be written in the 
frequency domain as x(a),z)—A( co)e The phase ve- 

locity 0 ) 1/3 is related to the phase (f> of the quotient of the 
Fourier transforms of the two time domain samples of the 
ultrasonic field at the axial locations Zj and z 2 according to 


c(io) 


0)(z 2 — Z 1 ) 
</>T 2N 7 7 


<f>(co) = tan 


Im(x( w,z, )/jt(o>,z 2 )) \ 
R t(x(u) 9 Zx)fx(<0 9 z 2 ))} 


(6b) 


With the geometrical constraints of the measurement system 
shown in Fig. 1, measurement of the ultrasonic field after 
propagation in the sample is only possible at z = L . The 
transducers placed in face-to-face contact, without the inter- 
vening wire specimen, is taken as the ultrasonic field at z = (). 
The integer N in Eq. (6a) is the integral number of wave- 
lengths contained in the length of the wire specimen, when </> 
is constrained to the range 0 to 277. 

From Eq. (6a), the phase velocities (c*j,c 2 ) and phases 
(<f)\,(f> 7 ) at closely spaced frequencies co { and <o 2 are related 
to N by 


CO 1 (f> 2 ~ CO 2 <f> | 

2 7r( (t) 2 CO] ) 

( C 2 C i ) [(f> ] <f>2 + 2N'7T((f> ] + (f> 2 ) + 4N~7T~] 

277 ( Z 2 Z i) (<O 2 -00\) 

(7) 


As discussed, the first longitudinal mode tends toward non- 
dispersiveness at low frequency. Nondispersive pulses propa- 
gate without distortion of the waveform envelope, and with- 
out change in the location of carrier frequency phase points 
within the envelope. A nondispersive or constant phase ve- 
locity approximation for N neglects the second right-side 
term. This approximation error term clearly decreases for 
smaller AT, which also occurs at low frequency. By transmit- 
ting a sufficiently low-frequency toneburst pulse, an accurate 
value for N at the burst center frequency can be determined 
using this approximation. 

Bursts are transmitted and the phase (f> measured while 
the center frequency of the pulse is incremented. Only the 
phase information at the dominant frequency of a particular 
pulse is used. By incrementing the pulse center frequency in 
sufficiently small steps, the phase wrap (jumps between up- 
per and lower limits of constrained phase) can be accurately 
tracked and the N incremented very precisely. Given an ac- 
curate value for N at the low starting frequency, N can be 
determined unambiguously over the entire frequency range 
of measurement. Thereby, the phase velocity is determined 
discretely as a function of frequency, using Eq. (6a). The 
range of measurement can be extended up to frequencies 
where interference from other excited modes occurs or where 
the pulse becomes so dispersive that interference between 
multiple reflections of the same mode occurs. 

This method was found to be more reliable than phase 
unwrapping, 25 wherein the constrained phase ((f)) is exam- 
ined from zero frequency to the frequency of interest. Each 
jump from 277 to 0 indicates that N be incremented by 1. 
However, only with the capability to move the measurement 
points for the phase spectral analysis into close proximity, is 
sufficient a priori knowledge available to determine N by 
phase unwrapping. When the measurement points are fixed 
and N necessarily is large, the density of the phase jumps 


256 J. Acoust. Soc. Am., Vol. 97, No. 1, January 1995 


Madaras et a/.: Dispersive pulse propagation in wire 256 



TABLE I. List of measured material properties. 


Wire type (diameter 
and density) 

£ ll (GPa) 

^LT 

G\ t (GPa) 

G tso (GPa) 

304 stainless steel 





406 fim (16 mil), 
7919 kg/m- 1 

211 

0.30 

74.5 

81 

305 fim (12 mil), 
7852 kg/m 1 
(annealed) 

208 

0.30 

76.0 

80 

127 /Am (5 mil), 
7913 kg/m 3 

225 

0.35 

71.1 

82 

aluminum 1100 





406 /urn (16 mil), 
2679 kg/m 3 

72.8 

0.36 

24.4 

27 

305 /xni (12 mil), 
2668 kg/m 3 

68.6 

0.36 

25.3 

25 

127 yum (5 mil), 
2742 kg/m 3 

71.0 

0.40 

26.0 

25 


increases. Most real dispersive pulses will have regions of 
frequency exhibiting nonphysically realizable phase (<f>) in- 
formation (i.e., corresponding to negative group velocity), 
often concomitant with low-magnitude complex amplitude 
A(o>). With a high density of phase jumps and even small 
regions of inaccurate phase information, the phase wraps be- 
come impossible to accurately track. 

As noted previously, curve fitting of measured phase ve- 
locities of the first longitudinal axisymmetric mode through 
variation of material properties allows complete material 
property characterization for waveguides of isotropic 
media. 22 The same method may be applied to transversely 
isotropic cylinders provided that E LL and v u are taken as the 
independent properties to be varied in obtaining the curve fit. 
Phase velocity measurements from low frequency to fre- 
quencies evidencing significant dispersion (approaching the 
cutoff of the second longitudinal mode) are required for ac- 
curate determination of both material properties. 

The lowest torsional axisymmetric mode is nondisper- 
sive, for both transversely isotropic and isotropic cylinders. 
This mode is sensitive only to G LT for transverse isotropy. 
Once the torsional wavepacket is identified, its velocity is 
measured by cross-correlation. 28 

Table I summarizes the properties obtained in this man- 
ner for six different wire types investigated. The final column 
for each wire gives the shear modulus which would be re- 
quired to support the conclusion of isotropy given the values 
of E ll and v LT determined by longitudinal mode measure- 
ment. Comparing this value of shear modulus with that ob- 
tained by torsional mode measurement is a measure of the 
material anisotropy and supports the conclusion that trans- 
verse isotropy better corresponds to the material behavior 
than isotropy, particularly for the steel wires. However, there 
is good agreement between solutions for the dispersion rela- 
tions using either Pochhammer’s exact solution for isotropic 
media taking E and v as independent properties, or using the 


numerical method described above and the complete set of 
five independent properties. The implication is that determi- 
nation of G lt was only of interest in verifying the anisotropy 
of the media, but this anisotropy has no significant effect on 
the first longitudinal mode behavior. Since the first longitu- 
dinal mode propagation was not sensitive to the exact values 
for EyY and these were set equal to E LL and <7 LT , 

respectively. 

Figures 4 and 5 show measured longitudinal mode ve- 
locities and the curve fits obtained, for 16-mil (406 /um) and 
5-mil (127 /urn) diameter stainless steel wires, respectively. 
The entire frequency range for which data was taken was not 
necessarily used for each wire. The solid circles indicate 
phase velocity measurements at discrete frequencies used in 
the curve fit shown as a solid line. In Fig. 5, additional mea- 
surements not used in the curve fit are shown as squares. 



Frequency (MHz) 

FIG. 4. Curve fit (solid line) of measured phase velocity (solid circles) to 
determine elastic properties for a 16-mil-diam 304 stainless steel wire. Ana- 
lytical curves due to 0.5% increase in E LL (dotted line) and 5% increase in 
v (dot-dash line) are also shown. 
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FIG. 5. Curve fit (solid line) of measured phase velocity (solid circles) to 
determine elastic properties for a 5-mil -diam 304 stainless steel wire. Ana- 
lytical curves due to 0.5% increase in E LL (dotted line) and 5% increase in 
v (dot-dash line) are also shown. 

Additional information regarding error analysis and accuracy 
is also conveyed and requires some explanation. 

The appearance of the measured dispersion at lower fre- 
quencies in Fig. 5 does not correspond to the first mode 
dispersion for any elastic cylinder. Such measurement 
anomaly was much more significant for the smaller diameter 
(127 /tm) wires tested. The frequency dependences of trans- 
mission through both of the transducer/wire interfaces, 
which has been neglected, may plausibly be influenced by 
the wire diameter. 

In general, several other factors may contribute to mea- 
surement error. First, the source waveform may be inad- 
equately characterized, due to the difference in source trans- 
ducer loadings between source measurement and 
transmission into the wire specimen. Furthermore, at the re- 
ceiving end of the wire, some of the wave incident on the 
transducer is reflected, which may also give rise to some 
phase shift. The frequency difference between lower- 
measurement frequencies and the resonant frequency of the 
transducer (5 MHz) may exacerbate these problems. 

If any resulting phase shift error at the pulse center fre- 
quency does not exceed some maximum value, say 7t/2, and 
the determination of N remains unambiguous, it is clear from 
Eq. (6a) that the phase velocity measurement becomes more 
accurate as frequency increases. The actual phase {<p) mea- 
surement becomes increasingly small relative to the quantity 
2ttN in the denominator of Eq. (6a). The error bars of Figs. 
4 and 5 indicate the change in measured phase velocity 
which would result from a ±7t/2 shift in (f>. 

Because of the apparent systematic error which occurred 
at lower frequencies for the 5-mil steel wire (Fig. 5), not all 
of the dispersive data available was used in the curve fit. 
That the analytical curve fit falls within the range of the error 
bars of the unused measurements for a tt/2 maximum error 
in <£, indicates that the phase shift error corresponds to less 
than a quarter of a wavelength. For the 16-mil diameter 
stainless steel wire (Fig. 4), measurements beyond 5 MHz 
were not used because of overlap and interference between 
longitudinal and flexural waves, as discussed previously. 

The dashed lines of Figs. 4 and 5 show the alteration of 
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FIG. 6. Pulse propagation of a 1-MHz toncburst in 25.27-cm-long, 16-mil- 
diam 304 stainless steel wire; (a) measured (solid line) and analytically 
extracted (dashed line) source waveforms, (b) measured (solid line) and 
analytically predicted (dashed line) receiver waveforms. 

the analytically determined curves due to 0.5% and 5% 
variations in E LL and i> LX , respectively. In comparing Figs. 4 
and 5, specifically the dashed and solid curves in each figure, 
it is apparent that in order to determine v u accurately, data 
associated with the increased dispersion at higher frequency, 
approaching the second longitudinal mode cutoff, must be 
available. Consequently, accuracy of the determination of 
v LT for the thicker wire was judged to be improved (i.e., the 
second mode cutoff is at much higher frequency for the thin- 
ner wire). The dispersion of the fundamental mode, as fre- 
quency and wire diameter to wavelength ratio increase, has 
been associated with the significance of radial inertia. 10,11 A 
conservative accuracy estimate of better than 0.5% for E LL 
and 5% for v LT is claimed when adequate dispersion mea- 
surements are available. 

VI. COMPARISON OF EXPERIMENTALLY MEASURED 
AND MODELED PULSE PROPAGATION 

Comparison of experimentally measured receiver wave- 
forms with analytically predicted responses, due to various 
measured sources, for the 16-mil steel wire, are shown in 
Figs. 6(b), 7(b), 8(b), and 9(b). The measured sources engen- 
dering these responses are shown in Figs. 6(a), 7(a), 8(a), and 
9(a) together with the source waveform extracted from the 
aforementioned measured receiver waveforms. The analyti- 
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FIG. 7. Pulse propagation of a 4-MHz toneburst in a 25.27-cm-long, 16- 
mil-diam 304 stainless steel wire; (a) measured (solid line) and analytically 
extracted (dashed line) source waveforms, (b) measured (solid line) and 
analytically predicted (dashed line) receiver waveforms. 

cally determined responses and sources have been rms nor- 
malized with the measurements they are compared with. The 
source excitations for Figs. 6, 7, and 8 would be described as 
relatively narrow band pulses with increasing center fre- 
quency (1,4, and 5 MHz, respectively), and with concomi- 
tantly increasing evidence of dispersion in the resulting re- 
sponse waveforms. These pulses are comparable to the 
pulses used to measure phase velocities and determine prop- 
erties. Correspondence in measured and analytical sources 
and receivers is reasonable, even for the highly dispersive 
propagation exhibited in Fig. 8. 

The efficacy of applying a semi-infinite waveguide 
model to propagation in samples which actually have bound- 
aries, and consequently reflected waves at the end z~L , 
should be explained. The measured response at z — L is due 
to both the incident and reflected waves, but only the inci- 
dent waves are being used in the propagation model. How- 
ever, provided that the reflected wave does not give rise to 
significant change in the phase of the measured response for 
a particular frequency harmonic and the reflection coefficient 
is constant as a function of frequency over the pulse band- 
width, the pulse waveform will not distort due to the pres- 
ence of the reflected waves. Similarly, the presence of mate- 
rial attenuation and the lack of consideration given to the 
material attenuation in the model does not give rise to sig- 
nificant discrepancies between the measured and modeled 


1.5 



-1.5 J 1 * 1 ■— 

47.0 52.0 57.0 62.0 67.0 

Time (microseconds) 



FIG. 8. Pulse propagation of a 5-MHz toneburst in a 25.27-cm-long, 16- 
mil-diam 304 stainless steel wire; (a) measured (solid line! and analytically 
extracted (dashed line) source waveforms, (b ) measured (top) and analyti- 
cally predicted (bottom) receiver waveforms. 

response, provided that the attenuation is relatively constant 
over the frequency range of the pulse bandwidth. Both the 
presence of reflected waves and material attenuation, under 
these conditions, affect the amplitude of the measured waves 
uniformly as a function of frequency. The appearance of the 
pulse waveform, within a multipicative constant, is not af- 
fected. Since the analytical response is rms normalized with 
the measured wave, these shortcomings in the model are not 
apparent. 

Figure 9 shows the pulse propagation due to a broad- 
band excitation of the source transducer resulting from input 
of an electrical spike. The pulse distortion is severe because 
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FIG. 9. Pulse propagation of a broadband pulse in a 25.27-cm-long, 16-mil- 
diam 304 stainless steel wire; (a) measured (solid line) and analytically 
extracted (dashed line) source waveforms, (b) measured (top) and analyti- 
cally predicted (bottom) receiver waveforms. 

of both the bandwidth and inclusion of higher-frequency 
components (transducer has 5-MHz resonant frequency). A 
consequence of this high level of dispersion is the lesser 
agreement between measured and predicted receiver wave- 
forms. Inexactness in the source measurement, used to pre- 
dict the response for a highly dispersive pulse, tends to be 
magnified in the response comparison. Reasonable agree- 
ment exists between measured and extracted source wave- 
forms, despite the lesser agreement in response waveforms. 
To some extent the same conclusion can be drawn from Fig. 
8, where the dispersion is lessened but still prevalent. 

Waveforms comparing measured and modeled propaga- 



(b) Time (microseconds) 


FIG. 10. Pulse propagation of a 5-MHz toneburst in a 25.22-cm-long, 5 mil 
diam 304 stainless steel wire; (a) measured (solid line) and analytically 
extracted (dashed line) source waveforms, (b) measured (solid line) and 
analytically predicted (dashed line) receiver waveforms. 

tion of a 5-MHz center frequency narrowband pulse, in the 5 
mil steel wire, are shown in Fig. 10. For the thinner wire and 
equivalent frequency excitations, much less dispersive pulse 
distortion is evident, reflecting the appearance of the disper- 
sion curves in Figs. 4 and 5. 


VII. CONCLUSIONS 

An analytical model of axially symmetric dispersive 
pulse propagation in cylinders with applied end conditions 
was experimentally investigated. Significant capability for 
modeling the generation and dispersive propagation of the 
fundamental longitudinal mode was demonstrated. Phase ve- 
locity measurement of the fundamental longitudinal and tor- 
sional axisymmetric modes were used to characterize mate- 
rial properties of the wires. Three of the five independent 
elastic properties for transversely isotropic media could be 
obtained, and the postulate of the nonisotropic nature of the 
drawn wires was supported. The method should be appli- 
cable to modeling the propagation of higher-order modes in 
solid cylinders, provided that the material is suitably charac- 
terized. 
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